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1 The functional determinant 



In order to get a feel for the spectral theory of natural differential operators on compact mani- 
folds, recall the idea of Fourier series, where one attempts to expand complex functions on the 
unit circle S 1 in C in the form 

oo oo 

Cfce lfce = Co + y^(aj cos jO + bj sin 

fc= — OO J = l 

The trigonometric series is an expansion in real eigenf unctions of the Laplacian A = —d 2 /d9 2 
(the eigenvalue being j 2 ). The exponential series is an expansion in eigenfunctions of the operator 
— \d/dft, which is a square root of the Laplacian; the eigenvalue is k. 

*This paper is a contribution to the Proceedings of the 2007 Midwest Geometry Conference in honor of 
Thomas P. Branson. The full collection is available at |http: / /www.emis.de/journals / SIGMA/MGC2007.html | 
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Suppose we wish to find the determinant of the Laplacian on the circle, or of its square root. 
Taking the naive view that the determinant should be the product of the eigenvalues in some 
spectral resolution, we immediately have some problems - we are really being asked to take the 
product of the nonnegative integers. We could declare this to be by convention (0 being one 
of the eigenvalues), but the question of computing a determinant like 

J ( d2 i 
det ("^ + 1 

still persists - is there any way to make sense of such determinants? 

There are, of course, other characterizations of the determinant of a finite-dimensional linear 
transformation T, and we might try to generalize one of these. For example, there is the exterior 
algebra definition, which appears in abstract index notation (see below for an explanation) as 

detT = T a \ ai T a \ 2 ---T a - am] . 

There is also a zeta function definition: let 

( T (s) = TrT- s , seC. 

(The exponential functions z z" s must be well-defined on spec(T); for this, since spec(T) is 
finite, we just need ^ spec(T).) The determinant is then 

detT = e'^. (1) 

Indeed, (JTJ) follows from the calculation 



d_ 

ds 



TrT~ s = -Tr log T. 

s=0 



The particular branch cut used to compute log T affects the value of Cr(0) (moving it by 
multiples of 2tt'\), but has no effect on its exponential. 

In fact, the zeta function approach is a fruitful one for operators like the Laplacian on S . 
Before pursuing it more deeply, however, let us look at still another approach. If T is symmetric 
and does not have in its spectrum, then 

detT = vr m f f e- {Tx ' x) dx] . (2) 



This is, in fact, the sort of formula that one tries to imitate in quantum theory, with the so-called 
functional integrals - integrals over function space. Since there are problems assigning meaning 
to such integrals, one tries to evaluate the infinite-dimensional analogue of the left side of (J2j) 
instead. Among the possible interpretations (or regularizations) of the determinant, one may 
choose the zeta function, and if the operator in question has certain properties, this succeeds in 
assigning a value. 

To see what zeta function regularization says about det(A) on S 1 , let us just brutally remove 
the eigenvalue. This makes the zeta function 

oo 

CAW = YJ 2 !' 23 = 2 CRiem(2s), 
3=1 

where CRiem is the Riemann zeta function. As a result, 
det A = e -<tic m (°) = (2vr) 2 . 
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If we could take the eigenvalue product interpretation of the determinant literally, this would 
say that the product of the positive integers is \/27r. 

oo 

Let us recall the process which assigns meaning to the expression CRi em (0)- The ser i es Yl 3~ s 

3=1 

converges uniformly and absolutely on the half-plane for Re(s) > 1. The resulting holomorphic 
function has an analytic continuation which may be constructed using Dirichlet series and the 
Mellin transform, to a meromorphic function on C, with a single simple pole at s = 1. 

What we would like to explore here is a generalization of this thinking to a situation in which: 

• S 1 is replaced by a compact n-dimensional smooth manifold; 

• the Laplacian is replaced by a differential operator with positive definite leading symbol. 

Roughly speaking, the same construction will go through in this setting, and there will be 
a functional determinant. 

For example, consider the Laplacian of an n-dimensional Riemannian manifold (M n ,g). One 
has the Riemannian connection V, which allows covariant differentiation of vector fields, and 
in fact of all kinds of tensors. If / is a smooth function, then VV/ is a 2-tensor. Contracting, 
one gets 

A/ := -g ab V a Vbf = -V a V a / 

in abstract index notation. Here indices which are repeated, once up and once down, denote 
contractions. Indices are raised and lowered using the metric, and the number and position of 
indices just indicate a tensor's rank and covariance type (rather than any choice of local frame). 
For more on abstract index notation, see |31| . 

The leading symbol of a differential operator D is, more or less, what one gets by writing D 
in local coordinates (and frames, if D is vector bundle valued), replacing each d/dx a by an i£ a , 
and keeping only the terms with the maximal number of £'s. Here £ is an indeterminate covector 
field, or one-form; the result is some sort of function on the cotangent bundle of M. A precise 
formula for the leading symbol may be given, without making choices of local frames, as follows. 
Suppose D has order k, and carries sections of a vector bundle V to sections of a vector bundle W. 
Then if / is a smooth function and y is a smooth section of V, 

a k (D)(x, {df) x )v = lim , 



t^oo t k 

where the subscript x denotes the value of a section (in this case, of the cotangent bundle or W) 
at x. Since any covector £ at x is (df) x for some /, this completely (and, as is easily checked, 
consistently) determines ak(D)(x,£)tp. 

For example, in the case of the Laplacian on ordinary functions, 

A/ = -g ab d a d b + (lower order), 

so the leading symbol is 

a 2 (A)(x,£)=g ab (x)Ub- 

The "2" in the o~2 notation just keeps track of the order of the operator (the maximal homogeneity 
in £). In this sense, the leading symbol of the Laplacian "is" the metric (or, more precisely, the 
metric inverse g ab ). 

Back in the general situation, o~k(D) may be thought of as a section of the vector bundle 



Hom(Symm(T*M <g> ■ ■ ■ ® T*M) <g> V, W), 
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using the identification of ^-homogeneous functions on a vector space (here, the cotangent 
space T*M) with symmetric fc-tensors on that space. The construction of a total symbol is 
another matter. One may view this as depending on choices (i.e., as a function carrying atlases 
of M which locally frame V and W to local total symbols) as in [23]. Or, one can assume more 
geometric structure and use it to try to define total symbols invariantly, as in [38] . This will not 
be an explicit issue here, though it is an issue in the underlying analysis of pseudo-differential 
operators, which provides the heat operators and complex powers of elliptic operators on which 
zeta function regularization relies. Roughly speaking, a total symbol will keep track of what 
happens to the operator under local Fourier transforms. 

The Laplacian's leading symbol has a precise positivity property: if £ 7^ 0, then 

g ab Ub > 0, 

by the positive definiteness of the Riemannian metric. We can speak of something like this 
in the vector bundle valued case, say for a A;— order differential operator carrying sections of 
a vector bundle V to sections of the same bundle V, replacing positivity by positive definiteness, 
at least if V has a positive definite metric. We say that D has positive definite leading symbol if 
for each £ 7^ in T*M, for each x G M, we have 

o~k{D)(x,£) positive definite in End(Vj;). 

A positive definite endomorphism A is, of course, one with (Av, v) > for each nonzero vector v, 
where (•, •) is the inner product (provided here by the metric on the bundle V). 

Note that only even-order operators can have positive definite leading symbol, since if k is 
odd and ak(D)(x,£) is positive definite, then af c (D)(x,—^) is negative definite. 

Perhaps the easiest examples of operators with positive definite leading symbol which do not 
have the form 

A h + (lower order) 

are the differential form operators Sd + add, where d and 5 are the exterior derivative and 
coderivative, and a is a positive constant not equal to 1. To be more specific, if <p is a p-form, 
then d(p is the (p + l)-form 

(dip)ao...ap = (p+ l)V[ ao <£ ao ...a s ...a p ], (3) 
and 6(p is the (p — l)-form 

(<fy>)aa...ap = -V ¥W..ap- 

In fact, one does not need a connection to define d; the right side of (|3j) is invariant as V ranges 
over all symmetric affine connections. 5 does depend on the connection V, and in fact is the 
formal adjoint of d in the form metrics 

This means that if ip is a p-form and rj is a (p + l)-form on a Riemannian manifold (M, g), then 
/ (d(p,rj) p+ idv g = / (<p,5r]) p dvg 

J M JM 

provided either ip or n has compact support. Here dv g is the Riemannian measure. The formal 
adjoint property makes it clear that 5 (unlike d) will generally change if we vary the metric g. 
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The form Laplacian 
Aform = Sd + d5 

is related to the Laplacian V*V (sometimes called the Bochner Laplacian in this context) by 
the Weitzenbock formula 

5d + dS = V*V + W. 

Here W, the Weitzenbock operator, is an action of the curvature tensor which has order as 
a differential operator. As a result, 

0"2(A form )0,£) = |£| 2 Idform, 

where the length is measured by the metric: |£| 2 = g ab CaCb- 
But if D has order k and E has order t, then 



In addition, formal adjoints (denoted by the superscript *) exist and are locally determined, and 



Thus 5d + ad5 has positive semidefinite leading symbol for a > 0, and for a > 0, the calculation 



(based on dd = 0) shows that 5d + ad5 has invertible leading symbol. Thus 5d + ad5 has positive 
definite leading symbol for a > 0. 

Beyond differential forms, in other natural bundles, the idea of positive definite leading 
symbol makes sense, but the calculus is not as simple - in particular, there is usually no complete 
analogue of the relation dd = 0. 

To construct the zeta function and the heat operator of a differential operator D on sections 
of V, we make the following assumptions. 

Analytic assumptions. (M n ,g) is an n-dimensional compact smooth manifold, and V 
is a vector bundle over M with a Hermitian metric h. A is a smooth-coefficient differential 
operator of positive order on sections V which is formally self- adjoint: 



and A has positive definite leading symbol. 

From these assumptions we may conclude that A has even order 21, since taking £ i— > — £ 
reverses the sign of the leading symbol for odd-order operators. In addition, L 2 (M, V) has 
a complete orthonormal set {(fj} of eigensections for A with real eigenvalues 



a k+i (DE)(x,0 = <T k (D){x,€)* e (E)(x,Z). 



a k (D)(x,0* =a k (D*)(x,0- 



(5d + add) (Sd + a^dS) = A| 



form 




A < Ai < • • • < Xj < ■ ■ ■ , 



satisfying 



A,- ~ const • j 



■2£/n 



J 



OO. 



(4) 



([!]) is the Weyl asymptotic law; see, e.g., [23 
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If r > 0, the Sobolev space L 2 (M, V) is the domain of the operator 

(A - A + l) r/2£ : ]T c m ^ Yl c ^ " A o + W 121 ^ 
j j 

in L 2 (M, V). The absence of any reference to A in the notation L 2 (M, V) suggests that this 
domain is independent of A, and in fact it is. (As a consequence of the compactness of M, the 
Sobolev spaces do not even depend on the metric of the tangent bundle or of V; see |28j.) In 
addition, there is always some usable operator A - for example, make metric and connection 
choices and take the Laplacian V*V. Though the Sobolev space L 2 (M, V) is Hilbertizable, it is 
more productive to think of it as a Banach space, with the equivalence class of norms represented 
by 

|M| 2 = f h((A-\o + l) r/2i ^(A-\ + iy/ 2e ^dv g . 
Jm 

This is just a more-or-less random choice of a norm from the relevant equivalence class. When 
delicate invariance considerations are in play, it sometimes pays to pick an invariant represen- 
tative of this class; see, for example, the last section below. 

The dual space is denoted L 2 _ r (M,Y); for real numbers r < s, 

(distributional sections) = £(M,V) C L 2 (M,Y) C L 2 (M,Y) 

C C°°(M,V) = (smooth sections), 

and L 2 = Lq. The Sobolev Lemma shows that 
C°° = f]L 2 r , 

r 

and this, together with the estimates defining the space of distributions, shows that 

r 

The heat operator e~ tA , t > 0, is of trace class, and is smoothing: it carries L 2 . continuously 
to l? s for any r,s, and thus carries distributional sections to to smooth sections. As a result, 
composition with e~ tA has a very "civilizing" effect on other operators: if, for some N, the 
operator B carries L 2 to L 2 _ N for each r, then Be~ tA and e~ tA B carry distributional sections 
to smooth sections. 

The heat trace has small-time asymptotics: 

Ti L 2e- tA ~ oo[A]far +ai[A]tsr- + ■ ■ ■ + a k [A]t^r- +■■■ , t [ 0. 

This means that if we take the sum on the right out to the k— term, we have accuracy to the 
order of the next term: 0{t^ n+k+l ^ 2t ). 

Note that one can conclude this accuracy immediately from a much weaker statement - 
accuracy to order 0{t N ^), where N{k) — > oo as k — > oo. To see this, pick k\ so that N(k\) > 
(— n + k + l)/2£, and expand out to order k±. Then the error in approximating Tv L 2e~ tA by 

— n —n-\-l —n-\-k 

aot~M + a\t~^~ + • • • + akt^r~ 

is 



t (-n+k+l)/2£ + . . . + akit (-n+ki)/2l + Q^.(-n+k+l)/2i^ = Q ^(-n+k+l) / 2ly 
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The way in which the ak depend on (the total symbol of) A can be minutely analyzed, and 
it is useful to do so. One qualitative observation is that the a Q( \d vanish in the current setting. 
(On manifolds with boundary, they consist entirely of boundary integrals.) Again, an excellent 
reference is [23J. 

In addition to the asymptotic series above, the heat operator trace also expands as a Dirichlet 
series of the first kind, 

oo 

Tr L2 (e- tA ) = J2 e ~ Xjt (*>0). 

3=0 

For example, when A is the operator —a?/d6 2 on the circle, the series is 



i + 5> 



3=1 

The corresponding Dirichlet series of the second kind take the form 



3=0 

Actually, this is not entirely accurate, since allows a finite number of Aj to be negative, and 
a finite number to be 0. Some sort of artificial "fix" is thus needed to write such a Dirichlet 
series; we shall make the choice 

A^O 

this is the zeta function Ca( s )- Whatever artificial convention we choose for dealing with non- 
positive eigenvalues, the effects will be felt in all succeeding formulas which make contact with 
the zeta function. 

The two types of Dirichlet series are related by the Mellin transform 



1 f°° 

(MF)(s) = — / t^F(t)dt, 
r (s) Jo 



under which 

M : e~ xt ^ A~ s (A > 0). (5) 
Because of flSJ), 

M : Tre^j, i-> T_£_^ 

= :Z A (t) ( A (s) 

when A has positive spectrum. When there are nonpositive eigenvalues, the convention on Ca(s) 
adopted above forces us to replace Za{£) by 

z A (t) = E e "* |Ajl = z A(t) ~ + 2 E ^h(\ 3 t), (6) 

Xj^O A J <0 

where q[A] is the multiplicity of as an eigenvalue of A. 

The analytic continuation of the zeta function, which so far is only well-defined on Re s > n/2l 
(by the Weyl asymptotics (J3J)), is accomplished by substituting the small-time heat asymptotics 
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into the Mellin transform expression for the zeta function. More precisely, assuming for the 
moment that there are no nonpositive eigenvalues, 



regular, Rcs>^ entim 



Since we can do this for arbitrarily large N, the IX function has been analytically continued to 
a meromorphic function with possible simple poles at the (n — k)/2£, where k is an even natural 
number. A useful viewpoint is that we are "approximating" the T£ function; the error is "small" 
when it is regular on a half plane starting far to the left. In this sense the "smallest" functions 
are the entire ones. 

To paraphrase the above, the T<^ function has at most simple poles on the real axis, at equal 
decrements starting with s = n/2£. Thus the zeta function has no poles at nonpositive integers - 
they're resolved by the zeros of 1/T(s). In particular, 

Ca(s) is regular at s = 0. 

This is what we needed in order to state the zeta function definition of the determinant, 

det A = e-&C°). 



Remark 1. The regularity of the zeta function at s = is a more delicate property than the 
Weyl asymptotics (|3|). It is possible to construct sequences of Xj satisfying the Weyl asymptotics 
for which the poles of the corresponding zeta functions perform various tricks - in particular, 
one can place a pole at s = 0. Somehow, regularity at s = indicates that an "organic" property 
of differential (or pseudo-differential) operators is being detected by the zeta function. 

Remark 2. If we don't necessarily have positive eigenvalues, the new spectral function Z^it) 
of (|6]) has its own modified heat coefficients dk, defined by 

oo 

-q[A] + 2^ a <o sinh(U,) + J>L4]^-"V 2 ~ £~ ^LAji^/ 2 . 

Note that d^ = for k < n. To build more determinant-like properties into the functional 
determinant, we might insert a sign to better monitor the negative eigenvalues: 

detA:= (_l)#{Ai<0} e -C^(0)_ 



2 Variation of the functional determinant 

We now begin to imagine wiggling the operator A in various ways. For example, we could let 
the operator A depend on the geometry (i.e. the metric g) like the Laplacian does, and vary the 
metric. An important way to vary the metric is conformally: 

g = n 2 g, 0<Oe C°°(M). 

In fact, by writing 

we can parameterize the conformal class of go hy to £ C°° (M) . 
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There are differential operators that are well adapted to this variation, namely the conformal 
covariants, or conformally invariant differential operators. Such an "operator" is really a rule 
assigning operators to metrics in a natural way, and has 

g = Q 2 g => Dtp = Q~ b D(Q a ip) (7) 

for some conformal biweight (a, b) (and all sections ip of the appropriate bundle). By rephrasing 
things in terms of operators on density bundles of appropriate weights, we could restate this as 
D' = D' , for a uniquely determined operator D', but we shall not pursue that viewpoint here. 
The best-known example of a conformal covariant is the conformal Laplacian 

n-2 

Y = A + — -K (K = scalar curvature), 

4(ra-l) 

on scalar functions; its conformal biweight is ((n — 2)/2, (n + 2)/2). 

The infinitesimal form of the conformal covariance relation is equivalent to the finite form ([7j , 
and is sometimes more useful. Let the metric run through the curve {g e u}e£R within a conformal 
class, and let 



d_ 

de 



£=0 



Then a conformal covariant has 

D'lp = aD(ivip) — bivDip, 

or, in operator terms, 

D* = -{b-a)uD + a[D,u\. 

(The final u in this formula should be interpreted as multiplication by u.) Infinitesimal confor- 
mal covariance is equivalent to finite conformal covariance, since any metric g^ in the conformal 
class may be connected to go by the curve {<?ew}ee[o,i] • 

In addition to the metric g, there are other standard ingredients involved in some geometric 
operators - the volume form E, and/or the fundamental tensor-spinor 7. If these are used, they 
should be scaled compatibly in forming D: 

e = n n E, 7 = srS. 

An example of a conformal covariant for which this is relevant is the Dirac operator y = 7 a V a on 
spinors. its conformal biweight is ((n — l)/2, (n + l)/2). Another example is the operator -kd on 
(n — l)/2-forms in odd dimensions; here the Hodge star operator carries p-forms to (n — p)-forms 
on an oriented manifold, and is defined by 

i*^) a p+ i . . .a n = ~^a x ...a n P 1 P ■ 

This makes E enter the definition of -kd on A^ n_1 ^ 2 ; the conformal biweight of this operator is 
(0,2). 

We would like to know what happens to our spectral data if we vary conformally, and if our 
operator A depends on g in a conformally "nice" way. It is generally hopeless to try to track 
the motion of individual eigenvalues (except in very special circumstances; see the discussion 
after (|10p below). However, it is sometimes possible to track the motion of quantities like the 
heat coefficients and the functional determinant. The reason is that these quantities have 
natural definitions in terms of traces of operators; such definitions take account of the whole 
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spectrum. If the variations of the operators are known, there is some chance of computing 
the variation of the natural spectral quantity. As we shall see, conformal variation of suitable 
operators yields a setting in which variational computations can be done. 
The good class of suitable operators is given by: 

Conformal assumptions. A = D h is a positive power of a positive-order conformal cova- 
riant D. 

For example, the conformal Laplacian satisfies both the analytic and conformal assumptions, 
since it is itself a conformal covariant. The Dirac operator y is conformally covariant, but 
does not have positive definite leading symbol. However its power y 2 , though not conformally 
covariant, has positive definite leading symbol (in fact, leading symbol |£| 2 Id), and so satisfies 
our analytic and conformal assumptions. The operator -kd on A^™ -1 )/ 2 does not have invertible 
leading symbol, so there is no hope of any of its powers satisfying the analytic assumptions. 

Under our analytic and conformal assumptions, we have the following: 

Theorem 1 (Theorem on Variations [151 EH USB [20] ) . With notation as explained just below, 
a k [A] m = (n- k)a k [u,A] = (n - k) / uU k [A]dv go , (8) 

J M 



CA0)' = ^j M U;^U n [A]-Y,\^ 



dv go , (9) 



where {fj} is any orthonormal basis of the 0-eigenspace of A. 



The [30. and |2Q] references really just treat the conformal Laplacian. 

The statement uses the local heat coefficients a k [uj, A] and C/fc[.A]. The first of these, a k [uj, A] 
is a term in the small-time asymptotics of Tr L 2Cue~ : 

j. a — n — raj-1 — n-\-k 

Tr L 2UJe tA ~ a [uJ,A\t— + a±[u, A)t^~ + ■ ■ ■ + a k [u;, A]t^~ + ■ ■ ■ , t[ 0. 

An analysis of the dependence of the heat coefficients a k [ui, A] on A and ui shows that they are 
integrals of differential polynomials in the total symbol of A, with coefficients that depend (in 
a universal way) smoothly on the leading symbol and linearly on the jets of a;. As a result, after 
integration by parts, 

a k [oj,A] = / ujU k [A]dv g 
Jm 

for some functions J7fcL4] built in a universal way from the total symbol of A. 
In particular, we have: 

Theorem 2 (Conformal Index Theorem [IS]). Under the above assumptions, the coefficient a n [A] 
is conformally invariant. 

Some other conformally invariant quantities which are important in these calculations are 

q[A] = dimAT(A) = #{Aj = 0}, ( A (0) = a n [A] - q[A], < 0}. (10) 

The invariance of q[A] is immediate from the conformal covariance relation. The invariance of 
#{Aj < 0} then follows from this and Browder's Theorem |18j . At the metric (fa) f° r some 5 > 0, 
the negative eigenvalues are contained in the interval (— oo, —5), and the positive ones in (<5, oo). 
With the invariance of the number of zero eigenvalues, this is good enough to keep the number 
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of negative eigenvalues constant on {g euJ }, for fixed uj and |e| bounded by some eq > 0. But this 
in turn implies constancy on the whole conformal class. 

The local data afc[u;,.A], C/fcL4] carry more information than afe[-A]; in fact 



a k [A]=a k [l,A] = J U k 



This more refined data (involving uj) is also spectral, but the spectra of many more operators 
(than just A) are now involved. 

If A is natural, that is, if it is built in a universal way from the metric (and possibly E 
and/or 7), then the total symbol of A must also be. But then, in turn, the £/fcL4] are also built 
in this way. By Weyl's invariant theory [37] . the £7jfcL4] are built polynomially from the Riemann 
curvature and its iterated covariant derivatives. 

For example, if K is the scalar curvature and A is the Laplacian on functions, 

U [A + aK] = (4^r n / 2 , 

U 2 [A + aK] = (^)- n / 2 {\-a)K, (11) 
U 4 [A + aK] = (4^)~ n / 2 • ^ • [90 (i - af K 2 - \r\ 2 + |i?| 2 - 30 (i - a) Ak\ , 

where R (resp. r) is the Riemann (resp. Ricci) tensor. One reason that formulas like this are 
important is hinted at by ©: we might be able to parlay explicit knowledge of C/ n [^4] into 
information about det A as a function on the conformal class {g u = e 2uJ go}. In fact, there are 
two genuinely different ways of doing this, one based on Q, and one based on both Q and ([S]). 

The first method is to simply integrate the variation along the curve {g eu) }, starting from go 
and ending at g^. As it turns out, the homogeneities are such that this is just the integration 
of a polynomial in e. Specifically, 

U k [A w ] = e~ 2k " (u k [A ] + E^AolM + ■■■ + [A ](uj)) , (12) 



(i) 

where Ut is a polynomial in the Riemann curvature and uj and their covariant derivatives, and 
is of homogeneous degree i in to. (In fact, it depends only on duo; see [6].) Here we have adopted 
the convention of labelling geometric objects depending on the metric by the subscript uj, rather 
than the subscript g u \ for example, 

Aui = A 9uj . 

On the other hand, 

dv u = e nuJ dv . 

As a result, U n [A £UJ ]dv £UJ is polynomial in e, of degree at most n. 

(By the conformal invariance of the conformal index jU n [A]dv, we may actually say "of 
degree at most n — 1". If n is odd, J7 n [^4] is identically zero, while if n is even, homogeneity 
considerations force the Ln n ' ) [A] term in (|12p to take the form c • |dw| n , for some universal 
constant c. But if c were nonzero, a simple scaling argument in uj shows that the conformal 
index is unbounded on the conformal class, contradicting the fact that it is constant.) 

The second method is an adaptation of the physicists' dimensional regularization. For simp- 
licity in this discussion, suppose that A has no zero modes (eigenvalues) at go (and thus on 
the whole conformal class). Suppose that we are interested in the functional determinant in 
dimension n. The idea is to get a formula for a n [A] which is valid in an infinite number of 
dimensions n' , including n' = n. Denote by Prim (for "primitive") the conformal anti- variation 
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corresponding to the conformal variation * above; this is well-defined up to a "constant of 
integration" - a constant function on the conformal class. The Conformal Index Theorem then 
says that 

f r , f + const, ri ^ n, 

Prim uU n [A] = { n ' n ^ 
J y 2i('a(®) + cons t> n' = n. 

The idea is now that the first formula continues to the second. To standardize the constants 
of integration, we view the primitives as functions of an indeterminate conformal factor u, 
measured from a base metric go, and require that our primitives vanish at go : 

/( a„[A^]-a n [A } / , 

1 ^(C^(O)-U(O)), n' = n. 

The analytic continuation to the special dimension n' = n is then expressed by the formula 

= ^(C^(o)-C^ (o)). 

n'=n 

As a result, we only need to know the integrated invariant a n L4], not the whole local inva- 
riant J7 n L4]. But we need to know it in "all" dimensions, or at least an infinite number of 
dimensions. If we have this information, we get a formula for the difference of log-determinants 
at two conformal metrics; that is, the quotient of determinants. 

Note that the first method (fixing the dimension and integrating along a curve in the con- 
formal class) is also giving us determinant quotients, since the definite integral described above 

computes C^,(0)-Cao(0). 

Of course, one cannot continue a formula from an infinite number of values of the independent 
variable to a special value. But such a continuation is possible in the case of a rational formula. 
That is, a rational function on C may be continued from any infinite set of values, because 
a polynomial function can. Alternatively, a meromorphic function on C for which oo is not 
an essential singularity may be continued from any set having oo (or any other point of the 
Riemann sphere, for that matter) as a limit point. What we need is a way of phrasing the above 
statements in terms of linear combinations, with rational coefficients, of a stable (under change 
of dimension) basis of local invariants. The ingredients out of which these invariants will be built 
will sometimes include the indeterminate conformal factor u - like the metric or the Riemann 
curvature, it will be thought of as just another formal variable, rather than a specific function 
on a specific manifold. 

Of course one knows that the heat invariants are usually not rational in the dimension. For 
example, for the operators A + aK on scalar functions of (jlip . it looks as if C/fcLA] will have 
the form (47r)~ n/2 V k [A\ with V k [A] rational in n (as long as the constant a = a(n) is rational 
in n). This is in fact true for all k, not just the k < 4 cases in (jlip . But the exponential 
function (47r) _n//2 certainly has an essential singularity at n = oo. The remedy for this is 
the fact that the (47r) -n / 2 factor is universal - one can factor it out of the coefficient of each 
invariant, and it's constant with respect to the relevant variation of the geometry. Thus it's 
really the coefficients of Vfe[.A] (and related quantities) that we're talking about, and these 
coefficients are rational. The possibility of factoring out a universal function of n so that we 
are left with rational functions of n is really a question about the form of the dependence of 
the operator's leading symbol on the metric. Similar results hold for other leading symbols; 
for example |£| 4 [10], and bundle-valued symbols which are not just powers of |£| 2 times the 
identity [3]. 



a n [A w ] - a n [Ao} 
n' — n 
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Demonstration of the Theorem on Variations. The key to the variational formulas is 
the conformal variation of the heat operator trace, or of the zeta function: 

Z{ty = (Tr e~ tA )' = -t Tr A'e~ tA = -t Tr D* hD^ 1 exp(-tD h ) 

aDui—buiD 

= th(b - a) Tr LuD h exp(-tD h ). 

Note that the step (Tve~ tA ) 9 = — tTr A'e~ tA is not valid on the operator level, even formally; 
the operators involved only have the same trace. There is no problem with domains due to the 
"civilizing" effect of the smoothing operator exp(— tA) - placing it next to a pseudo-differential 
operator of finite order produces an operator of order — oo, and a fortiori one that is bounded 
on L 2 . By homogeneity considerations, 

h(b-a) = 21, 



so 

(Tre- M )' = 2£tTTujAe- tA = -2it^- Tr ue~ tA . 
y ' dt 

However, 

+^L(+( k - n )/ 2£ \ — ^ ~ n Ak-n)/2i 

dt [t ] ~ 21 t 

So, provided that term-by-term variation of the asymptotic series is justified (see the remarks 
below) , 

a k '[A] = (n-k)a k [uj,A\. (13) 

From the zeta function point of view, the variational calculation is 

r(s)C(a)" = r(s)(Tr A~ s )' = T(s){Tv D~ hs ) 9 = I»Tr D 9 (-hs)D- fls - 1 
= (b- a)hsT(s)TrujD~ hs = 2£sT(s)C(lo, s). 

But a>k is the residue of the T( function at (n — k)/2£, so we recover (|13D . 
In addition, assuming for the moment that is not an eigenvalue of A, 



(s( A (u,s)) = 2£( A (u;,0). 



s=0 



But (a(u, 0) = a n \oj, A], since Res s= or(s) = 1. If is an eigenvalue for A, we carry out the last 
calculations after replacing exp(— tA) and j4~" s by the restrictions of exp(— t\A\) and |^4| _s to the 
orthogonal complement in L 2 of J\f(A), the null space of A. (Here a "restriction" is extended 
back to the whole space as on M(A).) On the level of kernel functions, we are then dealing 
with 

^ exp( — \\j\t)ipj (g) ip* and ^ \\j\~ s (fj (g) ip*. 

The zeta function arises from such a modified kernel, while the Z-function (as reflected in 
the C/fc[A]) arises from the "original" kernel function. Chasing through the consequences of this 
mismatch gives rise to the extra term in ([9]). 

In both the Z- function and zeta function arguments for the variational formula (JSj) -> there is 
an interchange of limit operations. In the case of the Z function, we are claiming that a certain 
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asymptotic series may be differentiated term-by-term with respect to an outside parameter (the 
e in the conformal factor e £u) ). That is, we are asserting that the asymptotics of the variation are 
the variation of the asymptotics. In the zeta function argument, an equivalent strong assertion 
is made: that the analytic continuation of the variation (again, with respect to the outside 
parameter e) is the variation of the analytic continuation. This is not just a matter of the 
uniqueness of the analytic continuation. In general, families of holomorphic functions with 
extremely nice dependence on an outside parameter may exhibit very ugly dependence in their 
analytic continuations. In either case, the issue comes down to suitable uniformity in the outside 
parameter of the estimates of error in the asymptotic expansion - in the zeta function argument, 
this enters when the heat operator and the Mellin transform are used to accomplish the analytic 
continuation. 

It's also worth noting that the issue of term-by-term differentiation of the asymptotic series 
with respect to an outside parameter (here, e) is more delicate than term-by-term differentiation 
with respect to the asymptotic parameter t. Indeed, since one can trivially integrate term by 
term, one can also differentiate: if 



/(t)~J> fc t fe , /'(t)~^M fc with 6 = 0, 



then 



f(t) + const ~£ J^+i, 



so that 



k + 1 



k + -1. 



This shows that the series for f'(t) is the term-by-term derivative of that for f(t), since the 
derivative of the term dot is 0. 

Note that the determinant quotient problem for conformal metrics is more or less trivial in 
odd dimensions (at least for manifolds without boundary), since ^7 n [A] = 0. (As noted earlier, 
the a Q dd vanish; the stronger statement that the U Q( ±d vanish is also true.) This means that the 
only variation in the determinant comes from the ip~ terms in ([9]). 

Here's a sample calculation to illustrate the dimensional regularization method of computing 
the determinant quotient. Our "input" operator A will be the square of the Dirac operator y = 
7 a V a on the bundle of spinors. (Recall the discussion just before the Theorem on Variations.) 
Suppose we would like the determinant quotient in dimension 2. Since 



a 2 [y 2 ] = -^2["/ 2 ](4^)-"/ 2 J K, 

the determinant quotient will be 
1 



n-2 



{{K dv) w - (K dv) ) 



(14) 



n=2 



If this is to work, there must somehow be a "hidden" n — 2 factor in the integral. To bring this 
out, recall the Yamabe equation 



n — 2 

Yu := Au + — —Ku 



n — 2 — n+2 



4(n - 1) 4(n - 1 

for the scalar curvature after a conformal change 

g = n 2 g, 0<neC°°, u = jrf. 



(15) 
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(Equation (|15p is actually the conformal covariance relation for the conformal Laplacian, applied 
to the function 1.) Rewrite (|15|) as 

. n-2 . n — 2 2=2, , 71 — 2 — n+2, , 

v ; 4(n-l) 4(n-l) 

. / "-a . , \ n — 2 n-2 , , n — 2 n+2 , , 

A e— ^ - 1 + —-Kqc~^~ = -K^e—". 

4(n — 1) 4(n — 1) 

On the second line, we have introduced a gratuitous —1 into the argument of Ao; this is harmless 
because AqI = 0. It is useful because it allows us to divide by (n — 2)/2: 

(n — 2 \ 1 n-2 , 1 n+2 , 

— -) + W^> K " e ^ = W^T) K ^ • (16) 

where 

e x - 1 



f(x) :-- 



x 

Evaluation at n = 2 now yields the Gauss curvature prescription equation 

A Lu + ^K = ^K UJ e 2 " (n = 2). 
In particular, 

{Kdv) = JiKdv)^ (n = 2), (17) 

since the difference of the two integrands, 2(Aouj)dvo, integrates to zero by the Divergence 
Theorem. 

(|17p just expresses the invariance of the conformal index in this setting. To get a handle 
on (|14p we need to go to higher order in the above calculation, in the dimension parameter 

0.^=1. 

2 

For this, first multiply (fl~6j) by dv$ = e~ ruj dv UJ : 
Now multiply by exp(/3u;): 
If we now set 

A" 



J 



2(n-l)' 
the above can be rewritten as 

(Jdv) u - (Jdv) = e^A {^—f^- ) + (e 2/3 " - 1)(J dv) 
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Now integrate and divide by n — 2: 

- y ( (Jd^ - (J A,) ) = J A j + Jo > d„ 



Evaluating at n = 2 (i.e., at /3 = 0), this gives 
1 



n - 2 



((J cfr;)^ - (J cfo) ) 



n=2 



-wA cj + U) J 



As a result, if the Dirac operator does not take the eigenvalue in the conformal class in which 
we are working, 



C^(0)-^,(0) = -^/{i.A„„ + u,J„ 



dvn 



Note that the objects that appear in the final result, A and J, really appear just because 
of the form of the integrated invariant a<i- In fact, since the only local scalar invariant that 
can contribute to a 2 is K (or J), the above formula is independent of A, except for the overall 
constant factor. For example, for the Laplacian on functions in dimension 2, 



Ci w (0) - Ca (0) = r^J { ^AoO, + u Jo 



dvn. 



We compute this by extending A for n = 2 to the conformal Laplacian A + (n — 2) J/2 in higher 
dimensions. 

3 Extremals of the functional determinant 

When we begin to think about maximizing or minimizing the determinant, the first thing we 
notice is that it is not invariant under uniform scaling of the metric: 

g = e 2a g, a£R A = e~ 2£a A, 

so 

Ca (s) = e 2£as CA(s), C' A (0) = &(0) + 2£aU0). 

Thus it's necessary to somehow penalize the determinant for the size of the metric. One way to 
do this is with a volume penalty: 

V(A,g) := (vol(7) 2£a(0)/n det A, 

is invariant under uniform dilation. (Recall that Ca(0) is a conformal invariant.) The formula 
for the conformal variation of the determinant says that the conformal variation of the 2?-func- 
tional is 

- log(V A )' = 2£vol(g)((uU n )~ - u>U n ), 
where the tilde denotes the average value. Thus the P-quotient is 
. V(A W ) 2/^(0). fe^dvo det(A,) 

log = log ^r- h log ■ 



V(A ) ' n 6 / du 6 det(Ao) 



Q-Curvature, Spectral Invariants, and Representation Theory 



17 



Now (to be somewhat vague) the coefficients in Ca(0) have something to do with the coeffi- 
cients in the determinant term - Ca(0) is a n LA], and the determinant term somehow came 
from a n [A]. As a result, though different A produce different P-functionals, combinations 
of a limited list of functionals keep appearing. For example in dimension 2, in the standard 
conformal class on the sphere S 2 , one is always looking at 

- log I e 2 ^^dt + / Lo{A u)d£, 
Js 2 Js 2 

where <i£ is the standard measure, normalized so that J g2 d£ = 1. This is exactly the quantity 
that appears in the celebrated Moser-Trudinger inequality. 
Recall the Sobolev embeddings 

L 2 ^ L 2n/(n-2r)_ 

The inequalities asserting the boundedness of the inclusion operators have sharp forms (more 
on this later; see ([28]) below) on S 2 . What happens as r — > n/2? One answer is the Sobolev 
lemma 

Lri. c — > . 
2 ~ l ~ t 

A deeper answer is the embedding 
T 2 ^ p l 

L n/2 ^ e • 

The inequality describing this inclusion compares an £ 2 / 2 norm with a functional that governs 
admission to the Orlicz class e L . When n = 2, this is the Moser-Trudinger inequality; for 
other n, it is the exponential class inequality of Beckner [5] and Carlen-Loss |19j . The sharp 
form of the Moser-Trudinger inequality is: 

Theorem 3 (Moser-Trudinger inequality), 
log f e 2 ^-^di < [ u){A u)d£, 

with equality iff e w is a constant multiple of a conformal factor on S 2 . 

To explain the term conformal factor as it's being used here, recall that there is a Mobius 
group of conformal transformations of the Riemann sphere. In terms of Riemannian geometry, 
such a transformation (diffeomorphism) has 

(h- l y 90 = n 2 h g . (is) 

The inequality is sharp because equality is attained for some u, and the case of equality is 
completely analyzed. 

In 1982, Onofri [27] used the Moser-Trudinger inequality to prove that on S 2 , the volume- 
penalized determinant of A is maximized exactly at the standard metric and its diffeomorphic 
images. Part of this is special to 2 dimensions - each metric on S 2 is diffeomorphic to one in 
the standard conformal class. Since diffeomorphisms preserve spectral invariants, it's therefore 
enough to study the standard conformal class. This is not the case in higher dimensions, even 
on the sphere. 

In fact, the space which may be loosely described by 

metrics 



Diffeo(M) k Conf(M) 
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is just a point when M = S 2 . (The semidirect product x in the denominator reflects the fact 
that diffeomorphisms act on conformal factors.) This space is generally a finite-dimensional 
object when n = 2 (in fact, Teichmiiller space), but it's a big, wild object when n > 2. 

One of the many questions that one can ask at this point is how all this plays out in higher 
dimensions - on S 4 , S 6 , and so on. Can we still solve a max/min problem by using sharp 
inequalties? What will the D-quotient functional and the sharp inequalities look like? In 
the following paragraphs, we look at the objects and methods that will generalize the two- 
dimensional versions appearing in the paragraph headings. 

The Laplacian. The Laplacian A on scalar functions turned up in the formula for the 
determinant quotient in dimension 2, even if it was some "exotic" A (like Dirac-squared) whose 
determinant we were measuring. In 1992, Graham, Jenne, Mason, and Sparling [25] proved 
the existence of a series of conformally covariant operators P m , of order m and of the form 
A m / 2 + (lower order); these operators exist for 

all even m when n odd, 
even m < n when n even. 

The Paneitz operator P4 was previously found in [29^ [32j [2Tj . 

In [39] , Wiinsch wrote a kind of formula for Pq - really more of an algorithm for computing 
it, but enough to prove its existence. 

The role of the 2-dimensional Laplacian in the determinant formulas is played by P n (the 
operator of dimension order), for even n. The role of the Gauss curvature is played by an 
invariant Q n (discovered in [16] in dimension 4, and in [8] in higher even dimensions) for which 

P m = P^ + 1 ^Y^ L Qm, P? n 1 = 0. (19) 

Dimensional regularization. The algebra of local invariants for general metrics can be 
simplified somewhat in the category of conformally flat metrics. These are characterized by 
the existence of local coordinates near each point in which the metric is g a t = Q 2 5 a b for some 
smooth positive function f2. There is an alternative characterization, in terms of invariants: if 
n > 4, a metric is conformally flat if and only if its Weyl tensor vanishes. The Weyl tensor, in 
turn, may be characterized as the totally trace-free part of the Riemann curvature tensor. One 
has the following formulas. Let 

j K 



2(n - 1) 
as above, and let 

n-2 

where r is the Ricci tensor. Then the Weyl tensor C is 

Cabcd = Rabcd + Hcffad _ Vbddac + V a d9bc ~ Vac9bd- 

As a consequence of the Bianchi identity, we have 

V a C abcd = (n - 3) (V c V& d - V d V bc ). (20) 

When n = 3, the vanishing of V c Vm — V^Vfe c characterizes conformal flatness (a fact which in 
itself is suggestive of some kind of dimensional regularization). In any case, for n > 3, in the 
conformally flat case, all local invariants are built from g, V and V (J being an abbreviation 



Q-Curvature, Spectral Invariants, and Representation Theory 



19 



for V a a ), and one has the "extra identity" that VV is totally symmetric. The ideal of identities 
(i.e., the ideal of left sides of identities with right side 0) is generated by ([2"0j) . together with the 
Ricci identities, which give formulas for [V Q , V&] applied to tensors of all ranks. 

We shall concentrate on dimensional regularization in the conformally flat category, though 
the restriction is not essential - there is a perfectly good algebraic setting in the category of 
arbitrary metrics. However, conformally flat metrics provide the cleanest route to the very clean 
extremal results on the spheres S even . 

Recall now our earlier formula for the determinant quotient: 

gca.w-^w)- '-'^:^ 

where n is the fixed dimension (in which we want a formula), and n' is the "running" dimension. 
What one gets out of this for even n > 2, given suitable dependence on n' , is a formula for the 
penalized determinant quotient functional of the form 

C (^77Tt[ (PnU,u)d£-log [ e n ^dc) + f {(B d£) u - (B d£) ) , (21) 

where B = B[A] is some local scalar invariant, and C is a constant. 

The invariant B is a linear combination of invariants from a precise (and remarkably short, for 
small n) list. For n = 2, this list is empty; for n = 4, there is one invariant that can contribute; 
for n = 6, there are 3, and for n = 8, there are 8. 

Remark 3. Besides being a convenient tool for calculation, dimensional regularization produces 
a theoretical result on the nature of the functional, which may be paraphrased as follows. If one 
follows the other, fixed-dimension method of integrating along a curve in the conformal class, 
then after subtracting out the correct multiple of the functional 

■FoM := ^-r I (Pnw,a;)de - log / e n ^d£, 
the dependence on u of the remainder is that of a local invariant: 

H") ■= [ ((B dOu - (B dOo) ■ 

We saw that in dimension 2, the penalized determinant quotient functional was an exactly 
the one estimated by the Moser-Trudinger inequality. What is the analogue for dimension 4, 
and higher dimensions? Remarkably, in the early 1990's, several lines of thought necessary 
for an understanding of this were maturing - higher dimensional Polyakov formulas [TB] , sharp 
inequalities [5j [19], and, though only an understanding of the Paneitz operator is needed in 
4 dimensions, the GJMS operators [25| . 

The sharp exponential class inequality in n dimensions is: 

Theorem 4 (Exponential class inequality). The quantity 



((P ft ) w,c;)d£-log / e<"-^dZ (22) 



2T(n) 

is > 0, with = iff e w is a conformal factor. 

The term conformal factor is again meant in the sense of (|18p . The analogue of the Mobius 
group, the conformal group of S n , is isomorphic to SO(n+ 1, 1); the action of this group will be 
needed in more detail later. 
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The operator (P n )o is the GJMS operator evaluated at the round metric; the precise form of 
this operator was given in [7J Remark 2.23]. Namely, let 



Then 

m/2 



n — 1 



( J P m )o = n{(^ + a-^ f B -a+^\\. (23) 

o=l v ^ , 

Recall that the GJMS operator implicitly defined a local invariant via (|19|) . Let 

n — m q n — m 

— 2 — Qm = PwXi so that P m = P m H - — Q m , 

where = 1. In general, P m and Qm are not unique - for example, one may add any multiple 
of the norm-squared of the Weyl tensor, |C| 2 = C abcd C a b c d, to the Paneitz operator P4 and 
still have a conformal covariant. What one really needs, for present purposes, is a choice of P m 
that is 

• has principal part A m / 2 ; 

• has conformal biweight ((n — m)/2, (n + m)/2); 

• has coefficients which are rational in the dimension n; 

• is formally self-adjoint; and 

• annihilates the function 1 in dimension n. 

Of course, such a P m is neither an operator on a particular manifold, nor an operator scheme 
for Riemannian manifolds in a certain dimension, but rather a formula, involving the dimension 
and some basis of the possible terms (built from g, V, and R) which is stable for large dimensions. 

There is certainly such a P m in the category of conformally flat metrics, and there is even 
an algorithm for writing it. Since such an operator must be A m / 2 for the flat metric g a b = S a b, 
we just need to write the operator invariantly for metrics (g^)a& = e2u ^ab- By the conformal 
covariance relation, 

(Pmh = e- (n+mW2 (P m ) e (n - mW2 . 

(As usual, we blur the distinction between a function / and the operator [if of multiplication 
by / - the final exponential in the foregoing formula should really be a multiplication operator.) 
In order to write this in terms of covariant derivatives in the metric g u , we write (P m )o as a com- 
position of conformal covariants evaluated at the metric go, and use the conformal covariance 
relations for these "smaller" operators. For example, 

(P m ) = A™ /2 = y m/2 = ( e (n+2)^/2 Fwe -(n-2)o;/2> ) m/2 ) 

where Y is the conformal Laplacian. As a result, 

= e -(n+ro)w/2/ e (n+2H2y we -(n-2) w /2\'n/2 e (T 1 -m)u J /2^ 

Another way to compute would be to use the fact that A = 5d, and write 
(P m )o = (M m/2 = (e^^e~("- 2 )-d) m/2 , 
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using the metric invariance of d, and the conformal biweight (n — 2, n) of 5. The result is 
(P m ) w = e -(n+m)w/2/ e no;^ e -(n-2)a;^'w/2 e (r l -m)w/2 > 

At any rate, we now commute all covariant derivatives to the right, past the expressions in uj. 
This produces terms in covariant derivatives of u. However, all second and higher derivatives 
of uj may be eliminated, in favor of local invariants, by iterated use of the Ricci tensor change 
formula 

1 

Uab = - Vab- U a U b + -UJ C UJ g ab 

and its iterated covariant derivatives. (Here ui a ...b abbreviates V& • • • V a a>.) The total degree of 
each term in e^ is 0, so we only have to worry about leftover occurrences of the first deriva- 
tive, Vw. There are none such, however, as a consequence of our having picked the correct 
principal part and biweight for the conformal covariant P m . The result is an expression for P m 
in the conformally flat category. 

In the case of high-order bundle- valued conformal covariants, a procedure like that just above 
is also possible, given the correct principal part and biweight. (The correct principal part may 
be deduced from spectral considerations [111 13].) In general, the principal part is written as a 
linear combination of compositions of "smaller" conformal covariants. 

The subject of promoting conformal covariants in the conformally flat category to conformal 
covariants for general metrics still features some open questions. It is known [24] that the 
operator on scalar functions with principal part A 3 does not promote in dimension 4; note that 
this is just beyond the restriction m > n on the GJMS operators in even dimensions. In terms of 
Bernstein-Bernstein-Gelfand resolutions (see, e.g., [22]), which classify (differential) conformal 
covariants in the conformally flat category, one knows that operators which are not longest 
arrows do in fact promote. Some longest arrows (like the critical order GJMS operator P n ) 
promote as well. A plausible conjecture is that the P n are the only longest arrows that promote. 

If we stay in the conformally flat category it is easily shown, by examining the above con- 
struction (which produces unique P m ), that the above conditions on P m hold. In addition [9], 
with the subscript still denoting evaluation at the round sphere metric, 

p / n+m \ 

(Qm)o = , n _ 2 m ^2\ ' and in particular, {Q n )o = r(n). 
1 \ 2 ) 

This means that the first functional in (|2ip . which we have called J-q{uj), is exactly the quantity 
estimated (i.e. asserted to be positive) by the exponential class inequality: 

W ■= TyfTT r / {P n u,u)d£ - log / e»( w - ffi >de. 
2T(n) J S n J S n 

The issue is now whether we can get an estimate of the other term, 
F{uj) := / ((BdOu-(BdOo), 

which is compatible with the exponential class estimate. This means that it should assert that 
CF{u) > 0, 

where C is the constant in (|2ip - the inequality should "go in the same direction as" the 
exponential class one. That is, it should assert this provided C / 0. If C vanishes, then Ct(0) 
vanishes; one way to see this is that the log term comes from the volume penalty; the case when 
no penalty is needed is exactly Gi(0) = 0. It is not clear whether this can happen in reality; it 
might be reasonable to conjecture that Cl(0) ^ for all A satisfying our assumptions. 
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Remark 4. It would also be nice if our compatible inequalities were compatible in another way - 
their cases of equality should contain the case of equality for the exponential class inequality. 
This will assure that the coupled functional vanishes exactly at conformal factors. But this 
additional property is provided automatically by the invariant nature of T{u) (Remark 3) - 
the functional is diffeomorphism- invariant, and thus in particular is invariant under conformal 
diffeomorphisms. This guarantees that when a; is a conformal factor, J-{oj) will vanish. 

To be more specific on how P n entered our formula for the J-q(u) part of the determinant 
quotient, recall first how the Gauss curvature prescription equation arose from the Yamabe 
equation. More generally, we would like a Q n prescription problem to arise in the same way 
from the Q m prescription problems, by dimensional regularization: 

Pm = P° m + t^Qm, O = 0, e^(P m U = (P m ) (e^f), 
n — m , ^ . n+m / n n — m „ \ n-m , , 



= (P° ) (e— w - l) + — ^— (Q m ) e— ' 
Dividing by (n — m)/2 and evaluating at n = m, we have 

(QnU™ = (P n °V + (Q m ) . 

In particular, since P® has the form 5Sd (by some invariant theory), 

{Qndv) u = \ (Q n dv)o. 



Going up an additional order in (3 = (n — m)/2, as in the m = 2 case, we get 
{Q n dv) u - (Q n dv)o 



= \J u(P n ) ujdvo + J uj(Q n ) dv . 

In this way, P n has naturally entered the determinant calculation. 
Here's the situation in dimension n = 4. P4 is the Paneitz operator. 

P 4 = A 2 + 5Td + ^^<34, 

where 

T = (n-2)J-4V; (V-tp) b = V a b <p a , Q 4 = \ J 2 - 2\V\ 2 + AJ. 

Since J C/4 [A] is a conformal invariant, L^LA], which a priori is a linear combination of K 2 , 
|r| 2 , |i?| 2 , and AK (or, equivalently, J 2 , | | 2 , |C| 2 , and AJ), must be something more special, 
namely a linear combination of Q4 , A J, and |C| 2 (C = Weyl). Since these quantities have 
conformally invariant integrals, the leftover functional 



/ 



((Bde) w -(Bde) )de 

can be expressed as a multiple of 

/ {{J*dO u - (J 2 ^)o) <ie. 

JS 4 
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(J 2 could actually be replaced in this formula by anything linearly independent of Q4 , A J, 
and \C\ 2 .) 

In higher dimensions, because of the Conformal Index Theorem, there is always some leeway 
in the choice in B. By p3], if L is a local invariant with j L conformally invariant, then in the 
conformally flat category, 

L = const • Pff + div, 

where "div" is an exact divergence, and "Pff" is the Pfaffian (Euler characteristic density). 
Since J Q n and J U n are conformally invariant, 

U n [A] = const • Q n + div, 

where the constant is uniquely determined by A. 

In dimension 4, the only possible universal divergence is AJ. Thus J 2 appears, basically, 
because it is a conformal primitive for AJ: 



■V 



JAlo 



(AJ)cj. 



As an overall conclusion in dimension 4, all of our penalized determinant quotient functionals 
end up being linear combinations of the two functionals 

JoM = \l wA (A + 2)u di - log [ e 4 ( w -*)de, 

3 J §4 J §4 

^(lo) = [ {(J 2 di) w - {J 2 dx) }. 

J § 4 

By the above, we know all about the functional Tq(uj). The other functional, J-\(u), is 
conveniently analyzed by the solution of the Yamabe problem [30l [36l E EH- The Yamabe 
quotient of the sphere §> 4 is 



y(f) 



[ /(A + 2)/de 



1/2 



The minimum value of y(f) is attained at / = 1, and so is 2. Thus for nowhere-vanishing /, 



2||/||1< jj\Ao + 2)fdt< 

by the Schwarz inequality, so that 
A / 



A / 



+ 2 



(24) 



2 < 



/ 



+ 2 



We would like to apply this to the nowhere-vanishing function 
This yields 



f = e 2 L 



4 < 



/ 4 ( e ~" ( ^i3 )e ") de 



This result is well worth stating in its own right: 
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Theorem 5. < / ((J 2 d£)w ~ (J 2 dx)o) , with = iff e^ is a conformal factor. 

The last part, about the case of equality, comes from the case of equality for the Yamabe 
problem on § 4 . The inequality really expresses the instance L\ <— > L 4 of the Sobolev Embedding 
Theorem in a way adapted to spectral invariant theory. As it has turned out, the functional 
J-\{oj) has exactly the same case of equality as the exponential class functional Tq{oj), not just 
a containing case of equality as described in Remark 4. 

Just what are the coupling constants between the functionals J-o(uj) and J-\(uj) for some 
prominent admissible operators A? Up to overall constant factors, here are the functionals 
when A is the conformal Laplacian, Dirac-squared, and Paneitz operators: 

Y: - IFoH - (w), 

y 2 : 33^o(^)+ 7^i (w), 

P 4 : 2l^oH-l&Fi(u;). 

Note that by dimensional regularization, these numbers really fall out of formulas for [A], 
provided we keep enough information to read the dimension dependence. (See |10j for the P4 
calculation.) 

Thus the extremal problem for T>(Y) and D(y 2 ) is solved, in the standard conformal class 
on § 4 . This calculation was originally done in [12], based on the determinant formula in [16], 

For the penalized det P4, functional, it's still unclear whether the round metrics are absolute 
extremals. The question can be phrased as comparing the sizes of the gaps (the differences 
between the left and right sides) in two sharp inequalities. The leading terms (those with the 
highest number of derivatives of u>) for Tq + aT\ produce a quadratic form that is 

positive definite for a < —8/15, 
indefinite for — 8/15 < a < —1/3, 

negative definite for a > —1/3. 

The number for P4 is —16/21 < —8/15; this is suggestive of a minimum for T>{(P^) U] ) at the 
round metrics, but the question is still open. 

A look at dimension 6 indicates that something interesting is going on. The P-functional in 
the standard conformal class on S 6 is a linear combination of four functionals, 

To(lo) = i / u;((P 6 ) u,)d£ - log / e 6 ^d£, ^{u) = [ (\dJ\ 2 d^, 
JF 2 (u;) = I ((|d/| 2 + 2J 3 )dO, - 54, 

= [ ((\dJ\ 2 + f J 3 - f J|F| 2 ) d£) w - 108. (25) 

Each may be estimated to be nonnegative, with equality iff e u is a conformal factor. Thus 
we "win", for a given operator A, if there are no conflicting signs in the expression of the 
P-functional for A as 

Remarkably, there are no conflicting signs in the cases of Y and y 2 : up to constant multiples, 
(00,01,02,03) is 

Y : (600,6,23,10), 

y 2 : -(11460,93,556,365). 

The nonnegativity of J-\{oS) is clear; it vanishes iff J is constant, iff (by Obata's Theorem) 
w is a conformal factor. An analysis of ^{oj) may be done by considering the Yamabe quotient 
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applied to the function J. The Yamabe quotient is 
/ u(A + 2J)u 



y(u) 



Hi 



so the Yamabe constant (minimum of the Yamabe quotient) on S 6 is 6 (setting J = 3 and 
u = 1). This gives 



({idj| 2 + 2J 3 }) w > 6 (| s6 (|j| 3 ( ia 



2/3 



since the Yamabe constant is conformally invariant (so in particular, it is the same at as 
at go). Now run an argument analogous to (|24p . based on (|28p below, and replacing the Schwarz 
inequality step with an application of the Holder inequality, with exponents | and 3. (A very 
general version of this argument is given in [9j Corollary 3.6].) The result is the nonnegativity 
of J~2(uj), with equality iff we are at a Yamabe metric, iff (by Obata's Theorem) u is a conformal 
factor. 

The estimate of J-^{uj) is more subtle. By [39J, there is a second-order conformal covariant D 
on trace-free symmetric 2-tensors. One may compute [9] that for n > 5, this operator has 
positive definite leading symbol, is formally self-adjoint, and has has positive spectrum on the 
round sphere. Since the sign of the bottom eigenvalue (if any) of a conformal covariant is 
a conformally invariant quantity, D has positive spectrum for metrics conformal to the round 
one. Thus for n > 5, 

(D UJ ip,p)(dv) UJ > 0, 

S" 

with equality iff (p = 0. If we substitute the trace-free Ricci tensor b^ of the metric g u for tp, 
the result (when n = 6) is the inequality J-^{oj) > 0. Furthermore, we have equality iff = 0, 
iff (S 6 ,^) is a space of constant curvature, iff a; is a conformal factor. 

An interesting aspect of the determinant calculation for the conformal Laplacian in dimension 
6 is the appearance of certain n — 8 factors. In fact, 

(4vr) n/2 7!a 6 [Y] = (n - 8) J {-3(n - 6)\dJ\ 2 - i(35n 2 - 266n + 456) J 3 
+ |(n - l)(7n - 30) J\V\ 2 - |(5n 2 -In- 48)trV 3 } dv 
= -§(n-8) J Q 6 dv + (n - 8)(n - 6) J {-f\dJ\ 2 

- ^(125n - 314) J 3 + |(7n - 5) J|y| 2 - §(5n + 28)trT/ 3 } du. 

This is analogous to an n — 6 factor that appeared in the 4-dimensional case. The root cause 
of this phenomenon is the relative conformal invariance of the heat invariant U n -2\Y]. To our 
knowledge, this was first remarked upon by Schimming [33], in connection with studies of Huy- 
gens' principle for hyperbolic equations. A Riemannian geometry proof is given in [30], based 
on the fact that U n -2 \Y] is the coefficient of the "first log term" in the radial expansion of the 
fundamental solution (Green's function) of Y. [T3] shows that this phenomenon is really a con- 
sequence of the conformal index property: take a fixed operator of the form D = — V a V a — E, 
and consider perturbations 

e - 2eu (D-6F), 
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where e and 5 are real parameters and u> and F are functions. Now compute 
d 2 _ 

' (*,e)=(0,0) 



d5de 



of the heat expansion with each of the two possible partial derivative orderings; the result is the 
relative conformal invariance of ?7 n _2[Y"]: 

U n . 2 [Y] UJ = e-^- 2 ^U n ^[Y] . 

(That is, the conformal variation of e^ n ~ 2 ^U n -2\y]w vanishes.) In the conformally flat category, 
this forces C/ n _2[^] to vanish. Thus, for example, £/g[Y] is forced to have an overall factor of 
n - 8. 

This calculation gives no special conclusion for powers of a conformal covariant (like y 2 ). 
Replacing Y by a higher-order operator leads to the relative conformal invariance of a different 
term in the heat asymptotics; for example, U n -^[P] in the case of the Paneitz operator P. 

Recently, Larry Peterson and the author [17] have computed the determinant quotient for 
the conformal Laplacian in the standard conformal class on S 8 . The starting point is a general 
formula for U$[D] due to Avramidi [2], where D is any operator of Laplace type. We use this to 
compute a$[Y] in arbitrary dimension n, in the conformally flat category. The result is (4-7r)~ n / 2 
times a linear combination (with rational-in-n coefficients) of the integrals of the nine scalar 
invariants 

(A J) 2 , J\dJ\ 2 , J|W| 2 , (V,dJ®dJ), J 4 , J 2 \V\ 2 , JtrT/ 3 , \V\ 4 , trT/ 4 . 

Here, for example, tr V 3 = V a bV b c V c a . Though this is a long calculation, something of a check 
is performed by the expected n — 10 factor - the fact that this factor appears is reassuring. We 
have not yet been able to use this to solve the extremal problem for T>(Y, g^), but several sharp 
inequalities are available. 

An interesting aspect of the extremal problem on low-dimensional spheres is the "checker- 
board pattern" 





dety 


dety 2 


s 2 


max 


min 


s 4 


min 


max 


s 6 


max 


min 



Here "max" means that the (suitably penalized) determinant quantity attains a max at the at 
the round metric and its conformal diffeomorphs. 

In (as yet inconclusive) work with Carlo Morpurgo and Bent 0rsted, we have attempted to 
explain at least the max/min alternation for the conformal Laplacian as follows. Suppose n is 
even. Morpurgo shows that for real s in the interval (5 



1, §), the zeta difference 



Cy u (s) -Cy (s) 



(26) 



is nonnegative. This difference has potential simple poles at the positive integers < n/2, so if 
all these poles are realized, and no zeros intervene, we have a quantity of sign (— 1) 2" 1 for small 
positive s; by the conformal index property, this should also be the sign of the log-determinant 
difference 



ck(o)-4(o). 



(27) 



Recalling the minus sign in the definition of the determinant as exp(— C'(0)), this agrees with 
the n = 2,4, 6 cases in the above table. 
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Considering what happens as we move along curves in the conformal class, by the obvious 
reflection principle, intervening zeros of (|26|) on the real axis must either come in pairs (i.e. be 
double zeros), or must move along the real axis. With double zeros, no harm is done to the 
above argument. A single zero, though, has the potential to reverse the sign of the above result, 
or ruin the result entirely, if it moves past s = 0. As it moves through s = 0, what we see 
is a zero of the log-determinant difference (|27h . So the problem is to rule these out - clearly 
a restatement of the assertion that the round metric is an absolute extremal. (Among the many 
things needed to make this discussion precise is to account for the conformal diffeomorphisms, 
which preserve the zeta function.) 

In any case, what is missing is some sort of general-dimension approach to the convexity 
properties of the determinant quotient. This would avoid the dimension-by-dimension approach, 
and have the side benefit of proving sharp harmonic-analytic inequalities rather than just using 
them. This is conceptually related to the topic of complementary series in representation theory, 
as we shall try to explain in the next section. 

4 The complementary series 

The constants and differential operators that appear in the sharp inequalities that estimate 
quantities like those in (|22p and (|25p . and in turn the determinant, also appear in the study of 
the complementary series of SOo(n + 1, 1). 

The conformal transformations h of any Riemannian manifold, ctran (M, g), form a group. 
Recall that such a transformation has 

h- g = (h- 1 y g = n 2 h g, o< n„ € C°°(M). 

The corresponding infinitesimal notion is that of conformal vector fields X; these satisfy 

C x g = 2u x g, u x £ C°°(M), 

where C is the Lie derivative, and form a Lie algebra cvf(M, g). The conformal factors satisfy 
cocycle conditions 

Qh 1 oh 2 = tth^hi ■ Clh 2 ), w [X!,x 2 ] — Xiujx 2 — XiVXi- 

The cocycle conditions are equivalent to the assertion that the family (series, in representation 
theoretic parlance) of maps 

u a (h) = n a h h-, U a (X) = C x + auj x (oeC) 

are homomorphisms into the group of automorphisms (resp. Lie algebra of endomorphisms) 
of functions on M, or in fact tensors, or tensor-spinors, of any type. The isometry subgroup 
(subalgebra) is defined by the condition = 1 (u>x = 0). 

When (M,g) is standard S n , these objects give the principal series of SOo(n + 1,1) (or 
Spino(n + 1, 1) if spinors are involved). That is, 

S n = G/MAN = K/M, 

where 

G = SO (n + l,l), M = SO(n), A = E+, N = R n , K = SO{n + l). 
The representations are Ind^ /j4Ar A ig) u 8) 1, with the interpretations 
A Rj tensor — spinor species; u w conformal weight. 
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The stereographic map M n S n is a special case of the more general Iwasawa decomposition 
theoretic map 

N ^ G = NMAN -> G/MAN. 

Such induced representations are best thought of algebraically at first, as (g, if)-modules. 
This means: take the iT-finite vectors, so one has a direct sum of K-modules. The action of q 
lets you travel among these modules - you can get from a if-module a to a .ff-module (3 only if 
5 ® a has a (3 summand in its if-decomposition (the selection rule). Here 

= t + S 

is the Cartan decomposition. 

An example is the expansion into spherical harmonics of functions on S n : 

£ — SO(n) Ej. 

Using the big Lie algebra so(n + 1, 1), we can move up and down the ladder, 

Ej-i <- -> 
within the spherical principal series representations 

ln d-M AN l ® ^® 1 =: (I7n +1/ ,u). 

The unitaries in this series are at imaginary v (the unitary principal series), and, on the real axis, 
in the interval \u\ < n/2 - the complementary series. The inner products on the complementary 
series representations (Un +U , u) have the form 

{(p,tp)u = / ip{A 2 v$)d£, 

where A 2v is the intertwining operator 

A 2v : (Un +1/ ,u) -> (Un_ V)U ). 

Since the two representations Un± v are naturally dual (they live in section spaces of the bundles 
of ^-densities), the integrand is natural. 

The Knapp-Stein intertwinor [26j is an integral operator, the integral converging in a certain 
range of Re(z/). Differential intertwinors, like the realization of the conformal Laplacian and 
Paneitz operator on the sphere, live in the analytic continuation (in v) of the family of Knapp- 
Stein intertwinors. In the sense of pseudo-differential operators and of complex powers of an 
elliptic operator [35j, assuming this analytic continuation, 

A 2v = A u + (lower order) . 

Up to a constant, (P m )o is A m . Not surprisingly, the obvious continuation of the formula (|23p 
given earlier for (P m )o turns out to be a formula for A 2u (all v): 

T{B + u + \) 

A 2 y 



where (recall) 



B = WA + 



n — 1 



2 
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For \v\ < n/2, by Stirling's formula for the Gamma function, the complementary series norm on 
U™ +u is a Sobolev l? v norm. 

The fact that the complementary series norm is an invariant Sobolev norm follows from the 
intertwining property: 

J ^ <f-A 2u { (c x + (| - f) oj x ) 4>}d£ = J ^ <p ■ (C x + (| + y) uj x ) A 2v i\) di 

Here we're using the fact that the case v = is in the unitary principal series; this really follows 
from the simpler fact that the volume form E has 

CxE = nuj x E. 

This last formula may be viewed as the origin of the rho-shift n/2. This also happens more 
generally in the theory of parabolic induction for semisimple Lie groups; the rho-shift due to 
the (g, 0) root system is related to the response of the natural volume form on G/MAN to 
infinitesimal transformations from g. In the present case, the calculation is 

0= f C x (<p-A 2u t/,dO= [J (£ x +(^ + v)ux)A 2u Adi 

*W l * : — ; — — — ' J 

+ J { (c x ~ Q + v) wx) MA ■ d£ + J (p(A 2l/ ip)nujx d£ 
= J ¥> j A 2 „ (c x + Q - v) ojx) Y> J<^ + J { (f-x + Q - vj lox) V?} A 2w i/j ■ d£. 

There's another natural invariant norm on the space carrying Ua_ VJ namely the LP norm 
with 

2n 

1 n + 2v 
Note that 

/-, \ ( n n\ 

P 6 (l,oo) ^ ^ (--,-). 

The invariance calculation is 

= fc x («) = P J <f?-\£ x <pW + / Vfmoxdt = pj ^ \C X + ^x) f 

this shows that when 

n n 

- = - + v, 

P 2 ' 

we have an invariant norm: 



/ 



The Sobolev embedding inequalities 

L^^L«-2^, L"+2" <^-> L_ v (v > 0), 



30 



T.P. Branson 



in the sharp form due to Beckner [5], compare these two invariant norms for v between —n/2 
and 0, saying that saying that 

max ; n (28) 

--v 

is attained exactly at conformal factors: both ip and ip should be constant multiples of £l£ , 
where h is a conformal transformation with (h~ 1 )*g = 0.\g. This gives a comparison of norms, 
and of norms on the corresponding dual spaces. These lead, by endpoint differentiation, to the 
exponential class inequality f)22[) . They are also used in deriving geometric inequalities like the 
estimate on ^(w) in ([25]) . based on [9j Corollary 3.6]. 

This interpretation in terms of representation theory is very clean and nice, but unfortunately, 
it has never been put to work in actually proving sharp inequalities. This remains a prospect 
for the future, much like (and probably closely related to) that of proving convexity properties 
of the determinant, and using them to solve extremal problems. The LP norms in (|28p are quite 
amenable to hard geometro-analytic methods like symmetric decreasing rearrangement, while 
representation theory is still quite linear (meaning quadratic in the present formulation - that 
is, concerned with inner products). 
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